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Applications of perturbative QCD to hard exclusive electroproduction processes in the Bjorken 
limit at small invariant momentum transfer t bring in a new type of parton distributions which 
have hybrid properties, resembling both the parton distribution functions/a ( x ) and the distribution 
amplitudes. Their i-dependence is analogous to that of hadronic form factors. We discuss general 
properties of these new parton distributions, their relation to usual parton densities and the evolution 
equations which they satisfy. 


INTRODUCTION 

The use of phenomenological functions accumulating information about nonperturbative long-distance dynamics is 
the standard feature of applications of pQCD to hard processes. The well-known examples are the parton distribution 
functions J p /h(x) Q used in pQCD approaches to hard inclusive processes and distribution amplitudes (p^fx) j^] 
etc. which naturally emerge in the asymptotic QCD analyses of hard exclusive processes. Recently, it was argued 
that the same gluon density f g ( x ) used for description of hard inclusive processes also determines the amplitudes of 
hard exclusive J/if> || and p-meson jj] electroproduction. It was also proposed || to use another exclusive process 
of deeply virtual Compton scattering y*p —> yp' (DVCS) for measuring parton distribution functions inaccessible 
in inclusive measurements. An important feature (noticed long ago |] ^|) is that kinematics of DVCS and hard 
elastic electroproduction processes, in a frame where the hadron moves fast, requires the presence of the longitudinal 
component in the momentum transfer r = p — p' from the initial hadron to the final: ry = £p. This means that 
the matrix element (p'| ... |p) is essentially asymmetric (“skewed”). For DVCS and p-electroproduction in the region 
Q 2 >> \t\,m 2 H , the “skewedness” parameter £ coincides with the Bjorken variable XBj = Q 2 /2(pg). For diffractive 
processes, both the skewedness £ of the nonperturbative matrix element (p'|...|p) and the absolute value of the 
momentum transfer t = (p 1 — p) 2 are small. Studying the DVCS process, one should be able to consider the whole 
region 0 < £ < 1 and t ~ 1 GeV 2 |J. In this situation, one deals with essentially nonforward (or off-forward in 
terminology of ref. jbj) matrix element (p'l-.Jp). The pQCD approaches incorporating asymmetric/off-forward 
parton distributions were formulated in refs, fqjj jl2| . A detailed analysis of pQCD factorization for hard meson 
electroproduction processes was given in ref. i| bib The basic idea of our papers 00 is that one should treat the 
initial momentum p and the longitudinal part of the momentum transfer r on equal footing by introducing double 
distributions F(x,y), which specify the fractions of p and ry carried by the active parton: k\\ = xp + pry. These 
distributions look like distribution functions with respect to x and like distribution amplitudes with respect to y. The 
longitudinal momentum transfer ry is proportional to p: ry = £p and it is convenient to parametrize matrix elements 
(p — r | ... |p) by nonforward TffX\t) or, in the t —» 0 limit, asymmetric distribution functions TffX) specifying the 
total light-cone fractions Xp = (a;-|-p£)p, (X — £)p s X'p of the initial hadron momentum p carried by the “outgoing” 
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and “returning” partonsj^. Double distributions F(x, y) are universal functions in the sense that they do not depend 
on the skewedness parameter £ while the asymmetric distribution functions J-q(X) form a family of X-dependent 
functions changing their shape when £ is changed^. The functions Fq(X) also have hybrid properties. When X > £, 
the returning parton carries a positive fraction (X — £)p of the initial momentum p, and F^(X) is similar to the usual 
parton distribution /(X). In the region 0 < X < £ the difference X — £ is negative, i.e., the second parton should 
be treated as propagating together with the first one. The partons in this case share the longitudinal momentum 
transfer r\\ = £p in fractions Y = X/£ and 1 — Y: in the region X < £ the function F^(X) looks like a distribution 
amplitude. 

It is very important to know spectral properties of the parton distributions F(x,y) and F^(X). In ref. Jl2| , it 
was shown that both arguments of a double distribution F(x,y) satisfy the natural “parton” constraints 0 < x < 1, 
0 < y < 1 for any Feynman diagram contributing to F[x,y). A less obvious restriction 0 < x + y < 1 @0] 
guarantees that the argument X = x + y( of the asymmetric distribution !Fq(X) also changes between 0 and 1. 
Since X = 0 can be obtained only if both x = 0 and y = 0, due to vanishing phase space for such a configuration 
the asymmetric distributions F^(X) vanish for X = 0. This property is very essential, because the hard subprocess 
amplitudes usually contain 1/X factors. When F^(0) ^ 0, one faces a non-integrable singularity Fq{X)/X at X = 0. 


ASYMMETRIC DISTRIBUTIONS 


There are two leading-twist quark operators z; A)?/> a (z) and ^ 0 ( 0 ) 757 ^( 0 , z; A)ij) a {z), where E( 0, z; A) 

is the usual path-ordered exponential which makes the operators gauge-invariant. In the forward case, the first one 
gives the spin-averaged densities f a {x) while the second one is related to the spin-dependent ones A f a (x). Here, we 
will discuss the ^ 7 ^ .5(0, z; A)ip a operators and gluonic operators with which they mix under evolution. The relevant 
nonforward matrix element can be written as 


{p', s' \ ip a (0)zE(0, z; A)tp a (z) \p, s)\ z 2 =0 (1) 

r 1 

= u(p', s')zu{p, s) J [e~ iX( - pz) F^{X-, t) - e ;(*-C) ( PUj?r?(X; dX 

+ u(p / , S l )^p u (p,s) £ (e~ iX ^!Cl{X-t) _ e <(*-0(i»)/c»(X;t)) dX, 

where z = z^y^. In Eq. quark and antiquark contributions are explicitly separated (cf. [p~5|). As emphasized by 
X. Ji ||, the parametrization of this nonforward matrix element must include both the nonflip functions Jy (X; t) and 
the spin-flip functions which we denote by /C£(X;f). Taking the 0{z ) term of the Taylor expansion gives the sum 
rules H 


[ |>c (X-,t) - X c a (X;t)] dX = F?(t ), 

Jo 

[ [Xl{X-t) - AC£(X;i)] dX = F£(t) 


( 2 ) 

(3) 


relating the nonforward distributions X“(X;t), /CjJ(X;i) to the a-flavor components of the Dirac and Pauli form 
factors, respectively. For gluons, the nonforward distribution can be defined through the matrix element 


<p' I (0)£“ 6 (0, z; A)G b a „(z) \p)\ z3=0 


= u{p’)zu{p) (z • p) 


1 

2 L 


D -iX(pz) 1 p i(X-Q(pz) 


F 9 c {X-t)dX+ ‘‘/Cf’ . 


(4) 

(5) 


In the formal t = 0 limit, the nonforward distributions X£(X';t), Af^(X;t) convert into the asymmetric distri¬ 
bution functions ■Ty(X), IC^(X). In the £ = 0 limit, F^(X) reduce to the usual parton densities X“ =0 (X) = 
/o(X) ; F 9 =q (X) = Xf g (X). For small £, T£{X) differs from / a (X) by 0(£). 


t Asymmetric distribution functions defined in ref. 0 are similar to, but not coinciding with the t —> 0 limit of the off-forward 
parton distributions introduced by X.Ji 0. 

*A more detailed discussion of double distributions and their relation to nonforward distributions was given in my talk about 

DVCS 10. 
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EVOLUTION EQUATIONS 


Near the light cone z 2 ~ 0, the bilocal operators .. ■ ip(z) develop logarithmic singularities In z 2 . and the 


limit 


0 is singular. Calculationally, these singularities manifest themselves as ultraviolet divergences for the 


light-cone operators. The divergences are removed by a subtraction prescription characterized by some scale y: 
!F${X\t) —* In QCD, the gluonic operator O g (uz,vz ) = z tl z v G < f a (uz)E ab (uz, vz\ A)G b av {vz) mixes under 

renormalization with the flavor-singlet quark operator. In the leading logarithm approximation, the easiest way to 
get the evolution equations for nonforward distributions is to use evolution equations [16 17| for light-ray operators: 




J-O a (0,z) = f [ B ab (u,v)O b (uz, vz)d(u + v < f)dudv, 
dy Jo Jo , 


( 6 ) 


where v = 1 — v and a,b = g, Q. For valence distributions, there is no mixing, and their evolution is generated by the 
QQ-kernel alone. Inserting Eq.fjcj) between chosen hadronic states and parametrizing matrix elements by appropriate 
distributions, one can get DGLAP and BL-type kernels and also to calculate the new kernels R ab (x, y; £, rf and 
Wf b (X, Z). The kernels R ab (x, y; £, y) govern the evolution of double distributions: 


d I" 1 r 1- ? 

y— F a ( x, y;t;y) = df V R ab (x , y; £, y) F b (£, y; f; y) dy . 

dy Jo Jo , 


(7) 


The form of equations is not affected by the t-dependence, and “t” will not be indicated in what follows. Since 
integration over y converts F a (x,y) into the parton density whose evolution is governed by the DGLAP 

equation, the kernel RQQ(x, y, £, 77 ) has the property 


fl-x 1 

/ RQQ(x,y;t,y)dy = - Pqq(x/£). 
Jo S 


Integrating Rqq(x, y; rf) over x one gets the BL-type kernel: 

fi-y 


Rqq( x > Vi (1 rf)dx = V QQ (y,rj). 


( 8 ) 


(9) 


Another set of kernels W£ b (X, Z) describes the evolution of nonforward distributions and asymmetric distribution 
functions: 


yjLf-iX; t;y) = f £ Wf b (X, Z) P b (Z-1; y) dZ . 
a V Jo . 


( 10 ) 


When £ = 1, the initial momentum coincides with the momentum transfer and J-^(X) reduces to a distribution 
amplitude whose evolution is governed by the BL-type kernels: WffL^X, Z) = V ab (X, Z). The BL-type kernels 
appear as a part of the asymmetric kernel W® b (X, Z) even in the general ( 7 ^ 1, 0 case [^2|. As explained earlier, if A' 
is in the region X < (, then the function P^(X) can be treated as a distribution amplitude H'c(T) with Y = X/f. For 
this reason, when both A' and Z are smaller than £, the kernels W“ b (X, Z) simply reduce to the BL-type evolution 
kernels V ah (X/f, Z/Q. Furthermore, the BL-type kernels also govern the evolution in the region corresponding to 
transitions from a fraction Z which is larger than ( to a fraction A which is smaller than (. Note, that the evolution 
jump through the critical fraction f is irreversible: X <( if Z < (. To put it in words, when the parton momentum 
degrades in the evolution process to values smaller than the momentum transfer fp = r, further evolution is like that 
for a distribution amplitude: the momentum can decrease or increase up to the r-value but cannot exceed this value. 

When X > £, the initial quark momentum Xp is larger than the momentum transfer r = Qp, and the asymmetric 
distribution function X“(X) is a generalization of the usual distribution function / a (X) for a skewed kinematics. 
Hence, the evolution in the region £<X<l,£<Z<lis similar to that generated by the DGLAP equation. The 
relevant kernels Pf b (X, Z) can be obtained from the light-ray kernels, e.g., 

P? Q (X, Z) = - £ B ab (w (1 - X/Z) ,w(l- X'/Z'))dw , (11) 

where Z' = Z — £. Note, that P ( ? < ^(A, Z) is given by a function symmetric with respect to the interchange of X, Z 
with X', Z'. In the f —► 0 limit, P®®(X, Z) converts into the DGLAP kernel: P^S 0 (X, Z) = pQQ (X /Z) / Z. 
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